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ABSTRACT 


As an important step in the characterization of a particular dynamic surface 
displacement transducer ( IQ I Model 501), this report presents one dimensional 
wave propagation in isotropic nonpiezoelectric and piezoelectric rods of variable 
cross section. With the use of the Wentzel-Kramers-Brillouin- Jeffreys ( WKBJ ) 
approximate solution technique, an approximate formula, which relates the ratio of 
the amplitudes of a propagating wave observed at any two locations along the rod 
to the ratio of the cross sectional radii at these respective locations, is derived. The 
domains of frequency for which the approximate solution is valid are discussed for 
piezoelectric and nonpiezoelectric materials. 
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INTRODUCTION 


In the characterization of a transducer designed for ultrasonic and acoustic 
emission measurements, it is important to understand the behavior of wave prop- 
agation in all components of the transducer, especially in the piezoelement. The 
piezoelement is important because it is through its behavior that ultrasonic waves 
are introduced into and received from the materials whose properties are sought 
by nondestructive evaluation ( NDE ). A transducers is usually designed with some 
specific measurement capabilities, such as 

(a) the type of waves (longitudinal, shear, surface, etc.) it can generate and/or 

detect; and 

(b) the frequencies at which it can best produce its desired response. 

In terms of its frequency characteristics, the transducer can often be classified as 
either narrow-band or broadband. The transducer whose characterization is sought 
in this study is a broadband type. The transducer is manufactured by Industrial 
Quality Incorporated (Gaithersburg, MD) and designated as IQI Model 501 . The 
sensing elements of the transducer consist of two components: a conically shaped 
piezoelement fabricated from lead zirconate-titanate ( PZT — 5A) and a backing 
material fabricated from brass. The transducer is designed for acoustic emission as 
well as ultrasonic signal generation. 

The aim of this report is to investigate how the piezoelement supports one 
dimensional longitudinal wave propagation. One of the earliest methods of ap- 
proaching this problem analytically was presented by Redwood [1]. Several other 
authors, including Filipczynski [2], Kazhis and Lukoshevichyus [3] and Williams 
and Doll [4], have subsequently dealt with the subject in various detail. However, 
unlike the present case where the piezoelement is conical, these authors dealt with 
cases where the piezoelement had a constant cross section. 

The first and the second sections of this analysis deal with the derivation of 
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the governing wave equation in isotropic nonpiezoelectric and piezoelectric rods 
of variable cross section, respectively. The third section is devoted to solving the 
derived wave equations by using a phase-integral method known as the WKBJ 
approximate solution technique. Through this method, an approximate formula, 
which relates the ratio of amplitudes of a propagating wave at any specific locations 
along the rod and the ratio of the cross sectional radii at these respective locations, 
is derived. The approximation is based on the first term of an asymptotic power 
series. Usually in such approximations, there is an inherent error that imposes a 
limit on the use of the approximate solution. These limits are also discussed in this 
report. 

Another problem with the approximation is the existence of points of discon- 
tinuity, known also as the transition points. Heading [5] summarized attempts by 
several authors to improve on the approximation so as to be valid at the transi- 
tion points. Gans [6], who studied the propagation of light in a stratified medium, 
examined the effects of the transition points in a systematic way. Since Gans, sev- 
eral authors have suggested ways to improve the WKBJ solutions. For instance, 
Bailey [7] produced an improvement when considering reflections of waves in an in- 
homogeneous medium, replacing the WKBJ solution by an approximation that is 
valid through a transition point. Bailey also discussed an iterative process whereby 
results could be further improved. Iami [8] produced a better representation of 
the WKBJ solution through the transition point and later applied the technique 
to heat transfer in a laminar boundary layer. Finally, Moriguchi [9] considered a 
phase-shift approach to a formula developed by Hines [10]. 

Concerning the validity of the WKBJ approximate solution, this report will not 
only address the conditions under which the approximation is true, but it will also 
establish the frequencies under which the conditions are met. Parameterized curves 
that relate the amplitude ratios to the frequencies of propagation are developed for 
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ONE DIMENSIONAL PLANE WAVE PROPAGATION IN 
NONPIEZOELECTRIC ROD OF VARIABLE CROSS SECTION 


In this section the one dimensional plane longitudinal wave equation is derived 
for an isotropic elastic rod of variable cross section. 

Consider an isotropic elastic rod of variable cross section extending from the 
origin z = 0 to z = L along the z-axis as shown in Fig. 1 . Assuming that a 
propagating plane wave becomes incident on the rod at the origin, an equation 
of motion which describes the propagating wave can be derived by considering a 
differential element dx of the rod as shown in Fig. 1 . 

A momentum equation of the differential element dx of the rod can be written 


as 


— aA(x) + ( a + 


E‘) 




= P 


A[x) + ~ — dx 
dx 


dx 

d 2 u(x,t) 

dt 2 


( 1 ) 


where A(x) [in 2 ] is the cross sectional area of the rod at any location along z-axis, 
p [lb-sec 2 /in 4 ] is the mass density, a [lb/in 2 ] is the longitudinal stress, and u(z,t) 
[in], which is a function of both space z and time t , is the displacement along the 
z-axis. 

Neglecting the terms of the order of dx 2 , eqn.(l) can be simplified to 




( 2 ) 


After rearranging the left-hand side of eqn.(2) as 


a 


dA(x) 

dx 


At \ da 

+ A(x) Tx 


dx 


[oA{x)\ 


( 3 ) 


eqn.(2) can be rewritten as 


M(*)] = pM*) 


d 2 u(x,t) 

~dt 2 


( 4 ) 
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The relationship between the stress a and the displacement variable u can be 


written via Hookes’s Law as 


a — E 


du(x, t ) 

dx 


where E [lb/in 2 ] is the modulus of elasticity of the rod. 
Substituting eqn.(5) into eqn.(4) gives 



A{x) 


du(x, t) 

dx 


pA(x) 


d 2 u(x, t ) 


dt 2 


(5) 


( 6 ) 


which is the one dimensional plane wave equation for an isotropic elastic rod of 
variable cross sectional area A(x ) . Note that E is a constant for homogeneous 
isotropic elastic materials. 

Eqn.(6) is a second order partial differential equation with a variable coefficient. 
An analytical solution can be obtained by a phase-integral method known as the 
WKBJ approximate solution technique. This technique is presented later in this 
report. 
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I 


ONE DIMENSIONAL PLANE WAVE PROPAGATION IN 
PIEZOELECTRIC ROD OF VARIABLE CROSS SECTION 


In this section, a one dimensional longitudinal wave equation is derived for an 
isotropic piezoelectric rod of variable cross section. 

The momentum equation (see eqn.(l)) is unaffected by the piezoelectric effect 
but the Hooke’s Law which requires an additional term [l] for piezoelectric materials 
becomes 


_ du(x,t) , _ 
a = E — t- 2 - 1 - hD x 
ox 


(7) 


Eqn.(7) shows that the presence of an electric flux D x in the z-direction pro- 
duces a “piezoelectric stress” D x , where D x [coul/in 2 ] is the electric flux density in 
the z-direction, and h [lb/coul] is the piezoelectric constant of the material measured 
with D x held constant. 

Substituting eqn.(7) into eqn.(4) and simplifying give 


E k 


A(z) 


du (x,t) 

dx 


-h 


dD x 

dx 


= pM x ) 


d 2 u(z,t) 

dx 


( 8 ) 


Even though the electric field is not uniformly distributed in the rod along 
the z-axis as shown in Fig. 2, there is no free charge inside the piezoelectric rod. 
Therefore, Gauss’s Law, which is given by [l] 


V • D = 0 , 


(9) 


where the vector D is the resultant electric flux density, remains valid. 
Eqn.(9) can also be written as 

dDx dD^y dDg _ 

dx dy dz ’ 


( 10 ) 


where D x , D y and D z repesent the z, y and z components of the electric flux 
density, respectively. 
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Since plane wave propagation is assumed, differentials with respect to y and z 
are zero. In particular, 


dD y dD z 
dy dz 


= 0 


( 11 ) 


Therefore, from eqn.(10) 



Considering eqn.(12), eqn.(8) reduces to 


( 12 ) 



du(x,t) 

dx 


pA{x) 


d*u(x,t) 

dt* 


(13) 


which is the one dimensional plane wave equation for a piezoelectric rod of variable 
cross section. 

Note that eqn.(13) is identical to eqn.(6) which holds for a nonpiezoelectric 
rod. Consequently, eqn.(l3) can also be solved by the WKBJ approximate solution 
technique which is presented in the next section. 
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W KBJ APPROXIMATE SOLUTION TO 
A SECOND-ORDER PARTIAL DIFFERENTIAL EQUATION WITH 

A VARIABLE COEFFICIENT 


The Method 

The method derives its name from the physicists Wentzel, Kramers, Bril- 
louin and Jeffreys, who were the first to develop and apply the technique to solve 
wave propagation equations for slowly-varying stratified media; hence, the name 
“WKBJ ” method. 

Stratified media are not relevant in the present case; however, the method 
will be applied to the wave equations derived for an isotropic rod of variable cross 
section. 


Application of the Method 

Recall the wave equation derived earlier for isotropic piezoelectric and non- 
piezoelectric rods of variable cross section: 


E l 


A{x) 


du(x,t) 

dx 


= pA{x) 


d 2 u(x, t) 

di 2 


(14) 


The WKBJ approximation technique presupposes a gradual variation of the 
cross section A(x) and assumes a solution of the form 


u{x,t) =a{x)e i9 ^ , 


(15) 


where a(x) is the amplitude of the propagating wave, i = \/^I , and 0 is the phase 
which can be defined as 

0(x,t) = f k{x)dx — u)t , (16) 

J *0 

where k(x) is the wave number of the propagating wave, oj is the circular frequency, 
and t is the time variable. 
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Eqn.(15) can de differentiated as 
du(x,t) _ d 


dx 


dx 


[ a(x)e' 9 


But a(x) is a function of x only. Therefore, 

SW*)I-5W*»I ' 

Considering eqn.(18), eqn.(17) can be expressed as 

= ^ {li + ' k 


du(x, t) 

dx 


Also, 

d 2 u(x,f) _ 5 
dx 2 

d 2 u(x,t) 
dx 2 

d 2 u(x,t) 
dx 2 

d 2 «QM) _ « 


= h [•" (£ + o(i) . 

= ik(x)c' e + tk(x)^j a(x) + 
= e'* [t'A:(x) ^ + ik(x) S j a{x) 


dx 


Tx + ** (l) 


ie a 

+ C dx 


Tx + ’* (x) 


a(x) 
a (x) 


= e " [(t + •'*(*>) (£ + “(*>) “I 1 ) • 


Similarly, 


~ ^ ~ = — P-^(x)w 2 a(x)e’^ . 


( 17 ) 


(18) 


(19) 


( 20 ) 


( 21 ) 


Substituting eqns.(20) and (21) into eqn.(14) and dividing through by e' 9 give 


dx 


+ ik(x) EA(x)^^ + ik(x)^a(x) = — pA(x)u; 2 a(x) . 


( 22 ) 


Solution of eqn.(22) can be obtained by introducing a dimensionless variable, 
£ , which can be defined as 


«-! ■ 


(23) 


where x is the independent variable along the direction of wave propagation and L 
is the scale length, which can be defined as the distance from the point of incidence 
of the propagating wave to the location of the output measurement. 
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Differentiating eqn.(23) gives 


— = I A 

dx L d£ 


( 24 ) 


Assuming that the wave number k 0 of the incident wave is large or simply that the 
frequency of the incident wave is high, the amplitude a(£) can be asymptotically 
expanded in powers of a dimensionless constant k 0 L by using Poincare’s definition 

in] 

«tf) = «o({) + j££«iU) + -+ + ■ ( 25 ) 

where ao(£) ><*i(£) > • • • ><*n(£) , • • • represent the amplitude terms of the asymptotic 
power series. 

Substituting eqn.(25) into eqn.(22) and rearranging give 


+ [pA(^)w 2 - £A(£)A: 2 (£)] |a 0 (£) + + ‘ * 


(26) 


= 0 . 

Multiplying the first, second, and third terms of eqn.(26) by [k 0 /k 0 ) 2 , (k 0 /k 0 ) 1 , 
and ( ko/ko)° , respectively, results in 


■ H + 

d£ k 0 L de 


■) 


'daotfl . 1 <*°i(0 , 

k 0 L d£ 


+ (^) l,to£ { A(f) * ({) l « 

+ T( [ EA({) * ({) (°° (f) + Ttfi + • • •)] } 

+ [pA(f)w 2 - £A(£)fc 2 (£)] ^Oq(0 + H 


( 27 ) 


= 0 
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Separating eqn.(27) by powers of 1/koL gives 
(j£j) = [M(e)w 2 -£:A(J)* 2 (£)]ao(e)=0 , 

(i^)' ; .*«{^(f)*(e)^ + ±[EA(()H()a, o(£)]} 

+ [M(Ow 2 -^(£)* 2 (f)]a 1 U)=0 , 

(£)’’ ^T ( [ EA ^ dJ W-} 

+ [ M (£) u . 2 - E J 4 ({)* 2 (£)] a 2 (£)=0 . 

The first expression of eqn.(28) gives the dispersion relation for the plane propa- 
gating wave, that is, 

M (£ V 2 -£ M (£)* 2 (£)=0 ( 29 ) 


or 


kHi) = Wo ■ 

But the phase velocity c of the propagating wave is defined as 


(30) 


c = 



(31) 


and so, considering eqn.(31), eqn.(30) becomes 

HO = 7 • (32) 

Since the circular frequency u and the phase velocity c are constants, the wave 
number k(£) is also constant and hence independent of £ . Denoting this constant 
as k yields 


k = - . (33) 

c 

If the dispersion relation as expressed in eqn.(29) is true, then the last terms of 
the second and the third expressions of eqn.(28) vanish. So, the second expression 
of eqn.(28) can be expanded as 

2 EMOHO^jf- +MOjf{EMOHO\=o ( 34 ) 
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and rearranged as 


[^K)W(f)] = o . 


(35) 


or 

£’A(^)A:ao 2 (^) = constant . (36) 

For any segment L of the rod between the x = x\ and x = x 2 , it can be 
deduced from eqn.(36) that 

a{x 2 ) 

_a(xi) 

where a(xi) and 0 ( 12 ) represent the amplitudes of the propagating wave at the loca- 
tions x = Xi and x = X 2 , respectively, and A(xi) and yl(x 2 ) are the corresponding 
cross sectional areas. 

Eqn.(37), which relates the variation of the amplitudes of the traveling wave 
to the cross sectional areas of the rod, was also derived by Green [12] and Liouville 
[13] who were the early investigators of this type of problem. Green [12] arrived at 
this solution by an energy conservation argument. 

Eqn.(37) can be simplified to 

0 ( 3 = 2 ) _ r(xQ 
a(xi) r(x 2 ) 

where r(xi) and r(x 2 ) are the radii of the rod at the locations x = xi and x = X 2 , 
respectively. 

Eqn.(38) is useful for quick practical approximations but it lacks information 
about the frequencies at which this approximation is true. In order to incorporate 
the information on the frequency into the approximation, it is necessary to inves- 
tigate the order of magnitude of the second term, ai(£) , of the asymptotic power 
series such that the approximate solution is true only when the first term ao(£) is 
retained. 



^(*1) 

A(x 2 ) 


(37) 
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The order of magnitude for ai(£) can be obtained from the third term of 
eqn.(28) which can be reproduced as 


k ° 2 r ( 


EA(() 


da o(£) 

di 


+ 




The second term of eqn.(39) can be expanded as 

<* o | [BA(£)a,(0|| 

= 2 £A(£)*^S + ai (()±[EA( £)*] 
Therefore, eqn.(39) becomes 

k ° 2 T( [ EA lf] + *> { 2EAk lt + Tt 1 ^* 1 } = 0 

where the argument £ is dropped for convenience. 

But from eqn.(34) 

£1**1 --(*54*)$ . 

Substituting eqn.(42) into eqn.(41) and rearranging give 


(40) 


(41) 


(42) 


‘4 


EA 


dap 

~dt\ 


+ 


2 iEAk(^--^^f) =0 

\ d£ op d£ ) 


(43) 


The expressions in the second set of parantheses in eqn.(43) can be rearranged as 

a 1 


da i ai da 0 d 

~dt ~ ~^~d£ = ~ a °di 


a 0 


(44) 


Substituting eqn.(44) into eqn.(43) and rearranging give 

£il = i (!n\( 1 ^ jL_ \EA da ° 

a 0 \ 2 \ k J \EAoq ) d£ [ 


d 

dt 


dt J 


(45) 


Since it is desired to have the results in terms of L , dividing both sides of 
eqn.(45) by k 0 L gives 

A r„_ 1 •• / 1 \ j r j- 1 

( 46 ) 


1 d 

ii 

_ i ( 1 > 

\ d 

FA da ° 

k 0 L d£ 

. 00 . 

2kL \ EAa 0 J 

* di 

L <K J 
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or 


' ai/koL' 

= *■ ( * ) 

\± 

\EA dao 1 

a 0 

2kL \ EAa 0 ) 

1 d( 

J 


( 47 ) 

Also since the only purpose here is to estimate how large ai/koL is with respect to 
ao , from eqn.(47) one can write 


d 

di 


OjJJcqL 

CLq 


=0 (*l) 


(48) 


where “O” represents “order of” . 

Note that since both sides of eqn.(47) have the same scale length, the term 


E Aclq d£ 

can be estimated to be order of unity. 
Integrating eqn.(48) results in 

a i/k 0 L 


EA 


dac 

dFJ 


(49) 


a 0 


-tt) ■ 


( 50 ) 


Substituting for £ from eqn.(23) into eqn.(50) and rearranging result in 

&-»(£)« ■ 


(51) 


Eqn.(51) gives the order of magnitude of the second term ai of the asymptotic 
power series as a function of the first term ao- 
Substituting for ai into eqn.(25) yields 


o(f) = «o(f) + O «o(f) + • • • + jj^«»(f) + 


(52) 


The concern here is to determine, analytically, the condition under which the 
higher order terms of the asymptotic power series, including the second term a \ , 
can be neglected in applying the WKBJ approximation. The condition can be 
derived from eqn.(52) as follows. Since a\ is of “order” ( x/L)/kL and also since ko 
was assumed to be large, then if 


x/I 

kL 


« 1 


(53) 
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the first and other higher terms of the power series can be neglected. Therefore, 
the WKBJ solution expressed in eqn.(38) would be true if the condition stated in 
eqn.(53) is satisfied. 

The physical signficance of the condition expressed in eqn.(53) can be obtained 
by expanding the equation as 

J2 « k ’ ( 54 ) 

or 

x 2tt 

L* K< Y ’ 

where A is the wavelength of the propagating wave. 

Consider a situation where the rod has a finite length L and where an observer 
is located at the end of the rod. In other words, the variable x is equal in magnitude 
to the scale length. In that case, eqn.(55) becomes 




(56) 


or 

A << 2 nL . (57) 

The inequality eqn.(57) is satisfied only by very short wavelengths or simply 
when the frequency of the propagating wave is high. Also, the inequality statement 
implies that the error associated with the WKBJ approximation is negligible at 
high frequencies. In short, according to Heading [5], the error associated with use 
of the WKBJ solution is expressed as 0(l/kL) times the approximation (eqn.(38)). 
Hence, 

■ M 

Eqn.(58) can also be derived from eqn.(52). Finally, the solution can be summarized 
as 


16 


A = R 


1 + ° [ kZ 


if 


7 «kL , (59) 

Li 

where A replaces the ratio a[x-i) / a{x\) y and R equals the ratio r(xi)/r(x 2 ) . 

The expressions of eqn.(59) represent the complete WKBJ solutions, after 
Wentzel [14], Kramers [15], Brillouin [16], and Jeffreys [17], for this problem. 

The conditional statement in the second expression of eqn.(59) indicates that 
the use of the formula is limited to cases when the frequency of the incident wave 
is high. How high the frequency must be is one of the issues of investigation in this 
report. The results are presented for some well-defined geometric parameters of the 
rod and are presented in the next section. 
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RESULTS AND DISCUSSIONS 


It is clear from the preceeding section that the validity of the WKBJ solu- 
tion rests on two important parameters, namely, the scale length, which in practice 
represents the length of the rod, and the frequency of the propagating wave. The 
WKBJ solution was derived under the assumption that the frequency of the prop- 
agating wave was high. How high the frequency must be in order to satisfy the 
solution is therefore an important question. 

In this section, investigation of the lower bound frequency at which the WKBJ 
solution is valid is pursued for both piezoelectric and nonpiezoelectric materials. 
One part of the investigation focuses on four piezoelectric materials and the other 
focuses on two nonpiezoelectric materials. 

The piezoelectric materials selected are lead mentaniobate ( PbNb^OQ ), lead 
zirconate-titanate (PZT — 5 A), quartz (St'02)> and lithium niobate (LiNtOz) with 
longitudinal wave speeds of [18] 131 x 10 3 in/sec (3,323 m/sec), 172 x 10 3 in/sec 
(4,350 m/sec), 226 x 10 3 in/sec (5,740 m/sec), and 291 x 10 3 in/sec (7,400 m/sec), 
respectively. These four materials represent the most widely used piezoelectric 
materials in ultrasonic and acoustic emission transducers. 

For this report, a range of values from 0.025 in (0.635 mm) to 0.25 in (6.35 
mm), which represents the usual thickness range of piezoelements of ultrasonic 
and acoustic emission transducers, was assigned to the length of the rod of each 
piezoelectric material. Six values, namely, 0.125,0.250,0.375,0.500,0.625, and 
0.750 , were selected as values of the ratio of the cross sectional radii. 

Based on the WKBJ solution 

A = R { 1 + u) ■ i <<kL ’ (60) 

where “O” is omitted for computational convenience, six plots (one for each value 
of R) were generated for each piezoelectric material. These plots are shown in 
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Figs. 3 through 6. Each plot represents a family of ten curves that relate the 
amplitude ratio A to the different frequencies of the propagating wave. Each curve 
was generated for one value of L . 

Before discussing the results of these plots, it is necessary to state that the 
WKBJ solution is satisfied when A is approximately equal to R , that is, 

A « R . (61) 

Although the frequencies at which the solution is satisfied can be deduced from 
Figs. 3 through 6 , it is probably not the easiest way to do so. For this reason, it is 
necessary to recall that the error associated with the solution (eqn.(60)) is given by 

Sm *{k) • < 62 > 

which in percentage terms may be written as 

£[%| = 100J2 (i) . (63) 

With this in mind, it becomes a matter of selecting an acceptable error and 
thereafter establishing the minimum frequency that satisfies the solution. In this 
report, am arbitrarily acceptable error of 1.4% is selected. The frequency at which 
the error is approximately 1.4% is termed the lower- bound frequency and is denoted 
by Ft . So, from eqn.(63), Ft can be defined as 

Ft ~ 16-— , (64) 

where c is the wave speed in the material, R is the ratio of the cross sectional radii, 
and L is the length of the rod. 

Eqn.(64) shows that Ft increases proportionally as R and decreases propor- 
tionally as L. Because R and L are variables, eqn.(64) is most easily understood if 
either R or L is a constant. 
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By fixing the value of L at 0.25 in and Ea at 1.4% , plots of Ft versus R are 
generated for the four piezoelectric materials and shown in Figs. 7 through 10 . The 
linear relationship between Ft and R can be observed from these plots. Similar 
plots can be generated for other values of L . 

In order to show the relationship between E and L, eqn.(64) can be reexpressed 
as 

E{%) « 16 R-^r . (65) 

rj* 

But the wavelenght Xt at the lower bound frequency can be expresssed as 

= jr ■ ( 66 ) 

Combining eqns.(65) and (66) gives 

E(%) « 16iE^ (67) 

So, based on eqn.(67) the relationship between E and L in parametrized form 
were generated and shown in Figs. 11 through 14 . For each piezoelectric material, 
the plots are generated for only 3 values of R because the relationship between E 
and R can easily be deduced from other values of R. It can be observed from Figs. 
11 through 14 that when the value of L is 0.25in, the value of E is equal to the 
acceptable error E a which was fixed to about 1.4%. A list of the estimated values 
of the lower bound frequencies and the corresponding acceptable errors are shown 
in Table 1 . 

Another aspect of this report is directed towards establishing similar param- 
eters in nonpiezoelectric materials. Aluminum and brass with longitudinal wave 
speeds of 250 x 10 3 in/sec (6,350 m/sec) and 175 x 10 3 in/sec (4,430 m/sec), re- 
spectively, were selected for the investigation. These two materials were selected 
primarily in preparation for experimental studies of wave propagation in isotropic 
elastic rods of variable cross section, which are presented in sbsequent report. Also, 
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brass was selected because its acoustic properties are very close to those of PZT—5A 
from which the piezoelement of the IQ I transducer was fabricated. For instance, 
the longitudinal wave speed in PZT — 5 A is only 1.7% lower than that of brass. 
Aluminum was selected because it is widely used in nondestructive testing ( NDT ) 
of materials and also for reference purposes. 

A range of values varying from 0.25 in to 2.5 in, which represents the usual 
thickness range encountered in UT of materials, was assigned to L . The same set 
of values of R was also selected in this case. 

With these parameters and using eqn.(60), parametrized curves that relate 
the amplitude ratio to the frequencies of the propagating wave are generated and 
shown in Figs. 15 and 16 for aluminum and brass, respectively. The interpretations 
of the results from these plots are quite similar to those presented above for the 
piezoelectric materials. 

Similarly, plots of Ft versus R are generated for aluminum and brass by using 

E a = 1.4% (68) 

in eqn.(64) but by assigning the value of 2.5 in to L . These plots are shown in Figs. 
17 and 18. Also, plots of E versus Xt/L are generated and shown in Figs. 19 and 
20 for aluminum and brass, respectively. The estimated lower bound frequencies 
and the corresponding acceptable errors are shown in Table 2 . Note that the 
lower bound frequencies are lower because the values of L for the nonpiezoelectric 
materials are ten the values of L for the piezoelectric materials. 

It can be observed that, in general, the lower bound frequency not only varies 
with respect to the geometric parameters of the rod, such as R and L , but also it 
varies with respect to the acoustic properties of the rod, such as wave speeds. This 
observation is important when applying the WKBJ approximate solution technique 
to solving wave equations for materials whose properties not only vary in cross 
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section but also in acoustic properties in the direction of the wave propagation, 
such as stratified media. 
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CONCLUSIONS AND RECOMMENDATIONS 


In this report, a formula that relates the ratio a(x 2 )/a(xi) of the amplitudes of 
a longitudinal propagating wave to the ratio r(xx)/r(x 2 ) of the cross sectional radii 
was derived by using the WKBJ approximate solution technique. Furthermore, 
the WKBJ solution was shown to hold only when the frequency of the propagating 
wave was high, that is, x/L << kL . 

How high the frequency should be was investigated for four piezoelectric mate- 
rials and two nonpiezoelectric materials. The minimum frequency, at which the ap- 
proximate solution is valid is within 1.4% error, was defined as the lower bound fre- 
quency. At the lower bound frequency, the error associated with the approximation 
was fixed at approximately 1.4% for all the materials discussed here. For instance, 
it was found that Ft was attained in all the piezoelectric materials when the ratio 
Xt/L was 1.4 , 2.8 , ... , 8.4 for the corresponding values 0.125 , 0.250 , . . . , 0.750 for 
R. So, for the piezoelement of the transducer with L = 0.100 in and R = 0.250, 
the lower bound frequency can be deduced to be approximately equal to 4600 kHz. 

In conclusion, one can say that the one dimensional wave propagation anal- 
ysis presented here not only enhances understanding of wave propagation in rods 
of variable cross section, but also forms a step in the characterization of conical 
transducers in general and, in particular, in the characterization of transducer IQI 
model 501 . 

Note, however, that the analysis neglects the effects of attenuation (energy loss) 
of the propagating wave in the materials and hence on the WKBJ solution. For 
this reason, it is recommended that the WKBJ solution presented here be verified 
experimentally. 
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Table 1 


Lower Bound Frequencies for Four Piezoelectric Materials 


Material 

Ratio 

R 

Lower bound 
frequency 

F T (kHz) 

A x 

T 

Error at 

m 


0.125 

750 

1.431 

1.400 


0.250 

1500 

2.862 

1.400 

Lead 

metaniobate 

0.375 

2250 

4.294 

1.400 

0.500 

3000 

5.925 

1.400 


0.625 

3750 

7.156 

1.400 


0.750 

4500 

8.588 

1.400 


0.125 

1000 

1.453 

1.400 

Lead zircon- 

0.250 

2000 

2.906 

1.400 

ate-titanate 

0.375 

3000 

4.359 

1.400 


0.500 

4000 

5.812 

1.400 

(PZT - 5A) 

0.625 

5000 

7.265 

1.400 


0.750 

6000 

8.718 

1.400 


0.125 

1300 

1.438 

1.400 


0.250 

2600 

2.876 

1.400 

Quartz 

0.375 

3900 

4.314 

1.400 

0.500 

5200 

5.752 

1.400 


0.625 

6500 

7.190 

1.400 


0.750 

7800 

8.626 

1.400 


0.125 

1600 

1.372 

1.400 


0.250 

3200 

2.744 

1.400 

Lithium 

0.375 

4800 

4.116 

1.400 

niobate 

0.500 

6400 

5.488 

1,400 

0.625 

8000 

6.868 

1.400 


0.750 

9600 

8.232 

1.400 


Fj 
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Table 2 


Lower Bound Frequencies for Two Nonpiezoelectric Materials 


1 

Material 

Ratio 

Lower bound 
frequency 

Ay 

Error at Ff 


R 

(kHz) 

L 

(%) 


0.125 

100 

1.389 

1.400 


0.250 

200 

2.778 

1.400 

Brass 


300 

4.167 

1.400 



400 

5.556 

1.400 


0.625 

500 

6.945 

1.400 


0.750 

600 

8.330 

1.400 

i 

! 

0.125 

143 

1.500 

1.400 


0.250 

286 

3.000 

1.400 

Aluminum 

0.375 

_ 

429 

4.500 

1.400 

0.500 

572 

6.000 

1.400 


0.625 

715 

7.500 

1.400 


0.750 

858 

9.000 

1.400 
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Amplitude ratio versus frequency for lead metaniobate rod with R = 0*375- 


R=0. 500, Material: lead metaniobate 



Fig. 3d Amplitude ratio versus frequency for lead metaniobate rod with R = 0.500, 


R*0.625» Material : lead metaniobate 
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Amplitude ratio versus frequency for lead metaniobate rod with R 35 0.625. 


R=0.750, Material: lead metaniobate 
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Ajnplitude ratio versus frequency for lead metaniobate rod with 
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R-0.375, Material: PZT-5A 
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PZT-5A rod with R = 0.500. 


R=0. 625 , Material: PZT-5A 
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Fig. 4f Amplitude ratio versus frequency for PZT-5A rod with R = 0.750. 
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Fig. 5c Amplitude ratio versus frequency for quartz rod with 


R-0.500, Mater ial : quart z 
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Fig. 5d Amplitude ratio versus frequency for quartz rod with R = 0.500. 
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Fig. 5e Amplitude ratio versus frequency for quartz rod with R = 0.625. 
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Fig. 5f Amplitude ratio versus frequency for quartz rod with R = 0.750. 
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Fig. 6b Amplitude ratio versus frequency for lithium niobate rod with R = 0.250. 




R=0.375, Material: lithium niobate 
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Fig. 6c Ampiiture ratio versus frequency for lithium niobate rod with R = 0.375 
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Fig. 6d Amplitude ratio versus frequency for lithium niobate rod with R = 0.500 



R=0. 625, Material: lithium niobate 
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Fig. 6e Amplitude ratio versus frequency for lithium niobate rod with 
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Fi^. 8 Lower bound frequency versus R for PZT-5A. 
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Kig. 9 Lower bound frequency versus R for quartz 



Material: lithium niobate 
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11a Error versus — for lead metaniobate rod with R = 0.125. 







-I 1 I I I I I Q 

OsJ O 00 C O CVJ O 

(%) 3 ‘J0JJ3 


59 


Fig. 12b Error versus for lead metaniobate rod with R = 0.250. 
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for lead metaniobate' rod with R = 0.375. 
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Fig. 12a Error versus — ^ for a PZT - 5A*s rod with R - 0.125. 
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Fig. 12b Error versus —2- for PZT - 5A rod with R = 0.250 
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Fig. 12c Krror versus — for PZT - 5A rod with R = 0.375- 
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Fig, 13a Error versus for quartz rod with R = 0.125. 
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0.125, Material: lithium niobate 
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Fig. 14a Error versus — for lithium niobate rod with R = 0.125. 
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Fig. 14b Error versus _ i- for lithium niobate rod with R = 0.250. 
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Fig 14c. Error versus — for lithium niobate rod with R = 0.375. 
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Fig. 15b Amplitude ratio versus frequency for brass rod with R = 0.250. 
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Amplitude ratio versus frequency for brass rod with R = 0.375. 
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Fig. 15d Amplitude ratio versus frequency for brass rod with R - 0.500. 
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Fig. 15f Amplitude ratio versus frequency for brass rod with R = 0.750. 
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Fig. 16e Amplitude ratio versus frequency for aluminum rod with R = 0.625. 
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Fig. 16f Amplitude ratio versus frequency for aluminum rod with 
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or versus — - for a brass rod with R = 0.375. 
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for aluminum rod with R = 0.250. 
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Fig. 20c Error versus 1 for aluminum rod with R - 0.375- 
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